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AN  EXACT  FORMULA  FOR  THE  MUTUAL  INDUCTANCE 
OF  COAXIAL  SOLENOIDS. 


By  Louis  Cohen. 


In  his  paper  on  Absolute  Standards  of  Inductance,  Dr.  J.  G. 
Coffin1  gave,  without  demonstration,  an  unpublished  formula  of 
Kirchhoff,  for  the  mutual  inductance  of  two  symmetrically  placed 
coaxial  solenoids,  with  the  remark  that  this  formula  is  an  exact  one 
and  very  valuable.  Although  the  calculation  of  inductance  by  this 
formula  involves  a  great  deal  of  labor,  yet  it  seemed  desirable  to 
calculate  by  it  the  mutual  inductances  of  some  solenoids  and  com- 
pare the  results  with  the  values  obtained  by  other  formulae  which 
are  only  approximate,  and  thus  get  an  estimate  of  the  degree  of 
approximation  of  the  other  formulae.  The  results  obtained,  how- 
ever, were  disappointing,  the  values  of  the  mutual  inductance  by 
this  formula  being  absurd  for  every  case  calculated,  thus  proving  the 
formula  to  be  seriously  in  error.  Feeling  quite  certain  that  Kirchhoff 
deduced  a  correct  formula,  which  was  probably  corrupted  in  passing 
through  several  hands,  Professor  Rosa  asked  me  to  undertake  to 
deduce  the  correct  expression,  and  so  ascertain,  if  possible,  just  where 
Kirchhoff's  formula  as  given  by  Coffin  is  wrong. 

Himstedt2  in  his  paper  on  the  determination  of  the  ohm  gave  an 
exact  formula  for  the  mutual  inductance  of  two  coaxial  coils, 
although  he  indicated  but  very  briefly  the  deduction  of  this  formula. 
It  occurred  to  me  that  the  formula  may  be  modified  so  as  to  give 
the  inductance  of  solenoids,  and  this  in  fact  proved  to  be  the  case. 
Before  proceeding,  however,  with  the  transformation  of  Himstedt's 
formula,  it  will  be  desirable  to  give  a  complete  deduction  of  the  for- 
mula, since  Himstedt  does  not  give  the  full  deduction  and  the  result 
is  of  considerable  importance. 

1  This  Bulletin,  2,  p.  125;  1906. 

2  Himstedt,  Pogg.  Ann.  28,  p.  338;  1886. 

295 


296  Bulletin  of  the  Bureau  of  Standards.  [Voi.3,No.2. 

The  mutual  inductance  between  two  circular  currents  is: 


M= 


II 


cos  e  dsdsl 


w 


■2li- 


rig.  1. 

If  the  two  circles  are  parallel  and  at  a  distance  d  from  each  other, 
then 

r*=A*-\-az+an—  zAa  cos  (0-00 

e   =  0—  0n  ds  =  ad<f>,  dsl  =  Ad(j)l 
A  and  «  are  the  radii  of  the  two  circles  and  therefore, 


M= 


Aa  cos  (0—0,)  d^dfa 


^Al-\-a%+dl—  2Aacos((f>—  0,) 
which  can  be  easily  put  in  the  following  form: 

^4<?  cos  ^rvA/r 


J/=47T 


J 


A/A*  +  a*  +  d*—2Aa"cosyjr 


If  2/  and  2/,  are  the  respective  lengths  of  the  two  coils,  d  the 
distance  between  the  centers  of  the  two  coils,  xx  and  x%  the  dis- 
tances of  the  centers  of  the  coils  from  the  origin,  then 


d=.\\  —  .\\. 
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and  the  mutual  inductance  of  the  two  coils  will  be  given  by  the 
following  expression: 


M=nn'  M  dx  dxx 


d-i     J -1      Jo  A 


47rnn'  Aa  C0S  ^d^  dx^dx*  (3) 

-y]A2  +  a2  -f  (xx  —  x2)2  -f  2Aa  cos  ^ 


n 


and  11'  are  the  number  of  turns  per  unit  length  of  the  two  coils. 
By  partial  integration  with  respect  to  ^r  equation  (3)  will  reduce 
to  the  following: 

A  a  sin  yjr  dx,dx„ 
M=A.irnri 


1 1 


-y/A2-\-a2-\-(x1— x2)2—  zAa  cos  -ty 


t    I  II  {Ad)2  sin2  yjr  d  yjr  dxtdx2 

+*irnn     I  I    [A*+a*  +  (x1-xi)*-2Aa  cos  ^]* 


(A a)2  sin2  i/raty  dx^dx2  , 

[A2-\-a2^r{x1-x,)2-2Aacos^Y  ^4) 


Integrating  twice  with  respect   to  xx   and  x2  and  substituting  the 
limits  we  get: 

M=^nn'     -X+s-^L  ^h?+(«H-'+4)'-  **»  cos  t 

Jo 

—  ^A2-\-a2-\-(d^-l1  —  I)2—  2Aa  cos  ^ 

—  ^M'+^+^+Z— /J2— 2^  cos  -^ 


+  A/^2  +  ^2  +  (^-/-/i)2-2^^  COS  i/r  (5) 

When  the  two  coils  are  concentric,  that  is  d=o,  equation  (5)  will 
reduce  to  the  following: 
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M=vcnn>    P,,  ,  2(faT^        \jA*  +  a*+t>- 2Aa  cos  ^ 
^  I   A%-\-ar—  2Aa  cos  -^[  v 

—  A/^"  +  ^2  +  ^i—  2^COStM  flf  t/t  T(6) 


where 


where 

V-- 


M=^mi\V-V1)  (7) 

J  77  

2(Aa)2sin2ylr^A2-{-a2-\-c2— 2 Aa  cos  ^  ,  , 
^2  +  tf2—2^  COS  i/r  ~~       ^ 


and  /^  is  obtained  from  /^  by  putting  cx  in  place  of  c.     The  evalua- 
tion of  Kcan  be  obtained  in  the  following  manner: 
Put 

</>=2^ 
then 

V 


S{Aaf 

J  2  sin2(/>cos2<£(y42  -f  « 2 + <r2  —  2Aa  -\-\Aa  sin2$) 

(^-f^2—  2^  +  4^*2  sin2</>)  ^A2+a2+c*—  2Aa  +  \Aa  sin^ 

Now  make  the  substitution  cos(£  =  jr  then, 

F=8(^tf)8  I  x2(i-x2){(A  +  ay+ci-4Aax*}dx 


dcf> 


Jo    I 


{(A  +  af-  4Aax2}^J(A  +  ^  )2 + ^2  -  4^^%/i  -  s? 


\dfw  C      x*-(i  +  k*)x*  +  k* 
+a)*\  (i  +  nxy(i-k*x*)(i 


8(Aa)*w  I         x*— (1 +  &)**  + JPx*      dx 

(A+a)*       (i  +  ^2y(i-W)(i-j2)  * 


where 


/&"=        AAa         n-  —  *~a    w*-(A4-aY4-c* 
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8(A 


I-\-7lX2—I 


+nx2)T](i-k*x*)(i-x*) 


1  —  n2x^  —  1 


afwVi    C 

K>2L7*  I  (1 

r 

J    (i  +  nx^ii-kVXi-x2)"' 

Ji  +  ^V-i  ,  "1 

(i  +  ;^2)V(i-^V)(i-^2)  ^J 

_8(^^)2^ji77_in      ;<i  +  /^2)  +  ^2   P 1 

-{A  +  df\n        n     +  7?  )j^k 


"*"    ;*2 


+  - 


<^jr 


—  nx 


k2x2)(i-x2) 


dx    (8) 


Now, 


«         «  1  , 


xidx 


1 


1  —  ;z.r 


^J{i-k2x2){i-x2) 


V(l-/^V)(i-^2) 


dx  =  F+7±E-—F 
k2         k2 


(9) 


i7,  iT,  II  are  the  elliptic  integrals  of  the  first,  second,  and  third  kind. 
It,  therefore,  remains  only  to  evaluate  the  integral 


j 


x*dx 


^j(i-k2x2)(i-x2) 


which  can  be  done  in  the  following-  manner: 


2^4 


^-*w-*}=^_*% 


2(\-\-k2)x 


.2 


V(l-^V)(I-^2)      '     -y/(l-k2X2)(l-X2} 

^k2xidx 


^dx 


r  2(i+#*)x*dx     c l 

I  V(i-^2)(i-^)+  I  J{i 


•.•^V(i-^v)(i-^) 

ldx 


-k2x2)(i-x2) 

dx 

-/&V)(i-P) 
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Hence 


; 


x'dx  -  -^  +  k*)E+  *±*F  (IO) 

V(i-AV)(i-*V  3*'  +   3*' 


Introducing  the  values  from  (9)  and  (10)  into  equation  (8),  and 
collecting  the  terms  containing  E,  E,  and  IT  separately,  we  shall 
finally  get 


V=\u 


(") 


The  elliptic  integral  of  the  third  kind  II  can  be  transformed  into 
the  complete  and  incomplete  elliptic  integrals  of  the  first  and  sec- 
ond kind,  by  the  following  relation:3 

k'  2  sin  6  cos,  6  I  J 

n  —  r\= 


V(i-sin26')(i-^sin26')| 

W+F{F(k'6)-E<ykf6)}-EF(k'6)  (12) 

where 

6  is  given  by  the  relation 

«=-i  +  £'2sin20  (a) 

F[k\6)  and  E(/c' \6)  are  the  incomplete  elliptic  integrals  of  modulus 
/-'  and  amplitude  6. 

From  equation  (a)  we  get 

u+i_{A*-ay  +  cXA-ay 
'  k'2  '~(A*-ay+c*(A  +  ay 

The  value  of  TI  as  obtained  from  equation  (12)  is  given  by  the  fol- 
lowing equation: 


SA.  Cayley:  Elementary  Treatise  on  Klliptic  Functions,  \  [83. 
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if  now  we  introduce  the  value  of  II  from   the  above  equations  into 
equation  (11)  and  simplifying,  we  shall  finally  get 

V=-(A*-a2)c[F\^F(k'0)-E(kfd)\-EF(k'd)\ 

S-(A2-6Aa  +  a*y-2(A*-a2yF 

3-y/(A  +  ay+? 

2(A*-\-a2)  —  c2  . _ _       ,  .2       2X7T  .     . 

V         3   '     ~4{A  +  df+<*E-c(A 2-a 2)-  (13) 

To  recapitulate  we  have: 

M=\Trnri(V-V?) 

V\s  given  by  equation  (13),  and  Vx  is  obtained  from  Vhy  replacing 
c  by  cx 

c—l-\-lx     ,     cx  =  l—lx 

F  and  E  are  the  complete  elliptic  integrals  of  modulus  k 

\Aa 


k2  = 


(A  +  af+c2 


F(kf)6)  and  E  (k' ,6)  are  the  incomplete  elliptic  integrals  of  first 
and  second  kind  of  modulus  k'  and  amplitude  6 

k^=T-k2 

■*«     (A2-ay+c*(A-ay 
sm  0- (A*__ay+(*(A  +  ay 

Formula  (13)  is  very  similar  in  form  to  that  given  in  Coffin's  paper. 
The  third  term  is  exactly  equal  to  the  corresponding  term  in  the 
KirchhofI  formula,  and  there  is  only  a  small  difference  in  the  first 
term.  There  is,  however,  a  very  great  difference  between  the  other 
two  terms  in  Kirchhoff  's  formula.  When  the  radii  are  equal  the  two 
formulae  will  become  identical,  because  the  terms  that  involve  the 
error  then  disappear.  The  test  that  Dr.  Coffin  applied  to  Kirch- 
hoff's  formula  by  deducing  from  it  the  self-inductance  of  a  solenoid 
by  making  the  radii  equal  is  thus  seen  to  be  insufficient. 
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The  accuracy  of  formula  (13)  is  shown  by  numerical  tests  given 
in  another  paper.4 

For  the  case  of  two  solenoids  of  the  following  dimensions: 

/=200,  /1  =  20.     A  =  i$1  a=io 

the  value  of  M  given  by  formula  (13)  and  confirmed  by  the  for- 
mulae of  Roiti  and  of  Searle  and  Airey  is 

M—\irnn'  X  6213.4. 

The  value  of  the  mutual  inductances  for  the  same  case  as  calcu- 
lated by  KirchhofPs  formula  as  given  in  Coffin's  paper  gives  the 

following: 

M—  —\irnn!  X  21399 

This  value  is  more  than  three  times  the  true  value  of  M.  If  we 
entirely  neglect  the  correction  factor  due  to  the  end  effects,  the 
mutual  inductance  is: 

M—A^a^^nn'  —  ^irnn'  X  6283.19, 

a  little  more  than  1  per  cent  greater  than  the  true  value  of  the 
mutual  inductance. 

The  value  of  the  self-inductance  of  solenoids  can  be  easily  obtained 
from  equation  (13),  as  was  done  by  Coffin,  by  making  A  =  a,  l—lv 
We  thus  obtain 

To  obtain  Vx  we  replace  c  by  cu  but  in  this  case 

cx  —  J—  /t  =  o 

*       3 

*This  Bulletin,  3. 
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Since,  however,  k?  =  r  A        ^ 5  =  1  in  this  case 

(A  +  ay  +  c* 


v  E—  I  2  ^1  —  sin2(£  dcj)=i 


v  F1  =  8^ 
1     3 


and 


,(^+4^V    _     4^2-^  .__ 2r^     8^31        ,     N 

z=4Hw?fe  +V-V4^+^"^)  (I4) 


Washington,  March  14.,  190J. 


